Route to turbulence in a trapped Bose-Einstein condensate 
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We have studied a Bose-Einstein condensate of 87 Kb atoms under an oscillatory excitation. For a 
fixed frequency of excitation, we have explored how the values of amplitude and time of excitation 
must be combined in order to produce quantum turbulence in the condensate. Depending on the 
combination of these parameters different behaviors are observed in the sample. For the lowest 
values of time and amplitude of excitation, we observe a bending of the main axis of the cloud. 
Increasing the amplitude of excitation we observe an increasing number of vortices. The vortex 
state can evolve into the turbulent regime if the parameters of excitation are driven up to a certain 
set of combinations. If the value of the parameters of these combinations is exceeded, all vorticity 
disappears and the condensate enters into a different regime which we have identified as the granular 
phase. Our results are summarized in a diagram of amplitude versus time of excitation in which the 
different structures can be identified. We also present numerical simulations of the Gross-Pitaevskii 
equation which support our observations. 

PACS numbers: 03.75.Kk, 47.27.Cn, 47.37.+q, 67.85.-d, 67.85.De, 67.85. Jk 



I. INTRODUCTION 

Turbulence is a phenomenon well known for classical 
fluids. Superfluids also exhibit this phenomenon [IHS] 
and in this case it is known as Quantum Turbulence 
(QT) . It consists of a collection of quantized vortices tan- 
gled in space [7]. In the case of superfluid Helium, QT 
has been achieved in many ways, for instance in thermal 
counterflow driven by heat injection [3]-[5] or by vibrating 
objects within the liquid [MTU]. 

It is believed that experiments with superfluids may 
overcome many of the complications encountered in clas- 
sical fluids [3HS] • In this sense, the existence of quantized 
vortices in quantum fluids may simplify the interpreta- 
tion and provide grounds for a theoretical description 
of the problem. However, due to the high density and 
strong interactions present in superfluid Helium, the size 
of the vortex cores is very small. As a consequence, the 
observation of QT and its dynamics is indirect and the 
interpretation of the experiments is difficult [8HT0] 

Trapped Bose-Einstein condensates represent an ideal 
system for studying QT. For being dilute and weakly 
interacting systems, the vortices can be directly observed 
by optical absorption. Additionally, because of the small 
scale and the high controllability that these systems offer, 
QT can be generated and studied under very controlled 



t Present address: 5. Physikalisches Institut, Universitat Stuttgart, 
PfafTenwaldring 57, 70550 Stuttgart, Germany 
'Electronic address: jseman@gmail.com 



conditions. Therefore, the emergence of QT in atomic 
superfluids may create new and exciting alternatives for 
exploring this phenomenon |llj . 

Recently, QT was detected in a BEC of 87 Rb atoms in a 
magnetic trap [12]. The turbulent state was achieved by 
introducing an external excitation through an oscillatory 
magnetic field. In a subsequent publication it has been 
discussed how the emergence of QT is related to the finite 
size of the atomic condensates [13] . The finiteness of 
the system is a particular property of ultracold trapped 
gases not present in superfluid helium. Therefore, even 
though QT is the same phenomenon in both BECs and 
superfluid Helium, the reason for its appearance in these 
systems can be very different. 

In the present letter we offer new studies in which we 
show how the variation of the parameters of the excita- 
tion must be performed in order to generate very different 
regimes in the sample. Specifically, we vary the time and 
amplitude of the excitation and observe the effect pro- 
duced in the sample. Depending on the combination of 
the time and amplitude, the excitation is able to gen- 
erate four different regimes, namely: (i) the bending of 
the main axis of the cloud, (ii) the nuclcation of regu- 
lar vortices, (in) the quantum turbulent regime, and (iv) 
the granular phase of the system, whose observation we 
report for the first time. Our results are summarized 
in a diagram of amplitude versus time of excitation in 
which the four different regimes can be identified. This 
diagram represents a very novel result because, due to 
the finiteness of the atomic superfluids, it is exclusive of 
BECs and, certainly, it is not present in superfluid He- 
lium. The diagram is also important because it shows 



the route in which the parameters of excitation must be 
combined to reach the different observed states. 

Of particular interest is the way in which the array of 
vortices evolves to the turbulent regime. For explaining 
this transition we have performed a numerical simula- 
tion of the Gross-Pitaevskii equation, obtaining qualita- 
tive agreement with the observations. Finally, we offer a 
plausible explanation for the granular phase and discuss 
future directions for its understanding. 



II. EXPERIMENTAL SYSTEM 

The experimental setup to produce the BEC and to 
excite it up to the turbulent regime is described in detail 
in our previous publications [12LI14J . In brief, we produce 
a cigar-shaped BEC of 87 i?6 with (2 - 3) x 10 5 atoms in 
a harmonic magnetic trap with frequencies uj x = 1~k x 
23 Hz and uj t = 2n x 210 Hz. Once the condensate is 
obtained an oscillating magnetic field is superimposed to 
the trapping potential. This oscillation is followed by an 
equilibrating time of 20 ms before releasing the atoms for 
free expansion. 

The oscillatory field has a fixed frequency of 2n x 
200 Hz and is produced by a pair of anti-Helmholtz coils. 
This excitation field generates distortions on the poten- 
tial, promoting a combination of deformation, displace- 
ment, and rotation in the cloud. As shown in references 
[12l - fT5] . the excitation is able to nucleate vortices in the 
sample and, for the proper conditions, to take it to the 
turbulent regime. 

The possibility of generating strongly noncquilibrium 
BECs with vortices and other coherent topological modes 
by modulating the trapping potential was discussed ear- 
lier [TB 18J. The prevailing creation of vortices by such a 
modulation owes to the fact that vortices with unit cir- 
culation are the most stable among all other topological 
modes. More precisely, the trap modulation, imposing 
no total circulation on the atomic cloud, creates vortex- 
antivortex pairs randomly distributed in space. This is 
different from the creation of vortices by an effective con- 
densate rotation 19, 20 . 



III. 



OSCILLATORY EXCITATIONS AND 
PRODUCED STATES 



We have considered a range of combinations of excita- 
tion times and amplitudes. For small amplitudes, inde- 
pendently of the excitation time, the only effect on the 
cloud is the bending of its main axis. That result was 
reported previously and it is an inherent consequence of 
the superfluid nature of the atomic cloud [22 152] . In- 
creasing the amplitude leads to a monotonic increase of 
the average number of formed vortices, although there is 
a big variation in the number of vortices when the same 
conditions are employed. The results for three different 
excitation times are presented in Fig. RTa). 
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FIG. 1: Average number of vortices observed in the excited 
cloud as a function of (a) the excitation amplitude for three 
different excitation times and (b) as a function of the excita- 
tion time for three different amplitudes of excitation. Lines 
are guides for eyes. The error bars show the standard devia- 
tion of the mean value of the number of vortices. 



The distribution of vortices in our sample does not cor- 
respond to any regular pattern similar to those reported 
in the studies related with the formation of vortex lat- 
tices [231 EI] . We believe that this is a consequence of 
our vortex formation mechanism, in which the nucleation 
of vortices and anti- vortices is equally likely. In fact, 
considering configurations of three vortices, we have col- 
lected evidence that it is vortices and anti-vortices that 
are formed in the cloud [5SJ . Even in the absence of reg- 
ularity, the states where the number of vortices in the 
sample can be well identified will be called regular. 

As the parameters of the excitation increase, when a 
number of vortices of the order of 20 is reached, a change 
in the behavior of the distribution of vortices in the cloud 
is observed. This corresponds to the formation of a ran- 
dom vortex tangle, one of the main features of QT. These 
distributions exhibit a large shot-to-shot variability. Ad- 
ditionally, the turbulent cloud demonstrates isotropic be- 
havior in its free expansion, as previously described [12] , 

The graph of Fig. [IV a) shows a clear connection be- 



tween the time and the amplitude of the excitation. Ex- 
citations for an extended period of time reach the turbu- 
lent condition for much lower amplitudes. Alternatively, 
the amplitude can be fixed and the time of excitation var- 
ied. Fig. [lib) shows the average number of vortices as a 
function of the excitation time for three fixed amplitudes. 
While no big difference is observed for the shorter exci- 
tation times, at longer times, the larger amplitude forces 
the system to the divergence in the number of vortices 
and, consequently, to the turbulent regime. 

After reaching QT, increasing even further the time 
and/or the amplitude of excitation renders the sample to 
a complete granulation, where the pieces of condensate 
are spread out over the cloud. The grains persist if the 
pumping of energy is kept flowing into the sample. The 
granular state of BEC is a kind of a heterophase mixture 
where BEC droplets are surrounded by normal gas |26j . 
Experimental confirmation concerning the superfluidity 
of each grain still remains to be addressed. 

The diagram of Fig. [2] shows the correspondence be- 
tween the observed distributions and the combinations 
of time and amplitude of excitation. It exhibits four do- 
mains, each corresponding to one of the observed regimes. 
The border lines between these domains are just guides 
to eyes. This diagram is useful for understanding the 
route to the turbulent regime. It shows that the im- 
portant quantity related to this route is the product of 
the amplitude and excitation time. This is equivalent to 
considering the total energy pumped into the cloud and 
could be used for characterizing the threshold behavior. 
In fact, the number of formed vortices is related to the 
total energy pumped into the system. The threshold for 
reaching turbulence also must depend on the size of the 
cloud, as discussed in Ref. [13] . Since QT occurs when 
the number of vortices is around 20, one should expect 
that the cloud must be sufficiently large to contain many 
vortices. 



IV. NUMERICAL CALCULATIONS FOR 
TURBULENT REGIME 

For better understanding the experimental results, 
we have performed numerical simulations based on the 
Gross-Pitaevskii equation, with the addition of a phe- 
nomenological term, 7 to account for dissipation. In 
dimensionless units, the equation for our simulations is 
given by [271 EH] 
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The net potential acting on the atoms is the sum of the 
harmonic magnetic trap and the oscillatory field, and can 
be approximately expressed by 
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FIG. 2: Diagram showing the domains of the parameters as- 
sociated with the observed regimes of the atomic superfluid. 
The points are the experimental observations. Figures on the 
top correspond to typical observations. For the region of regu- 
lar vortices, the number of vortices varies with the parameters 
as presented in Figs. [If a) and nib). Gray lines are guides for 
eyes, separating the domains of different observations. 

+ [y cos $i — x sin 9\ — 82 (1 — cos f2 £)] 
+ [z cos 62 + x sin 02 — 83 (1 — cos Qot)] > , 

where A = u x /ui r and t = A4 (1 — cos$7 £) are time de- 
pendent angles. For our experimental conditions, f^o = 
2ttx200 Hz, A x ~ tt/60 and A 2 ~ tt/120. The amplitudes 
for the translational oscillation of the potential minimum 
are {61,62,63) = a(2, 5, 3)/zm/a r , where a r = yjh/mu) r 
and a is a variable parameter that represents the am- 
plitude of the center-of-mass oscillation. Hence, the pa- 
rameter a is proportional to the amplitude of the excita- 
tion. We employ the representation 'F = ^p{y, z)cf>(x) in 
Eq. (Til) and consider 2D simulations in the y—z space. In 
this case, the interaction term becomes U2D = 4araN/R x , 
with R x being the characteristic size of the condensate 
along a;— axis. For our experiment U2D = 740. Since 
the thermal atoms, which cause the dissipation, move 
together with the potential, we also have to consider the 
reference frame co-moving with the potential. In this 
frame Eq. (fTl) becomes 
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with the linear momentum p = — zV and the angular 
momentum L = f x p. We consider Vb 



(y 2 + z 2 )/2, 



v(t) = (O,v y ,0)smn t and fl(t) = (fl x , 0, 0) sinfV- Us- 
ing the half of the oscillation period T — tt/Hq, we ob- 
tain v y ~ 2S2/T = 2^82/^ ■ The rotation frequency Cl x 
is also estimated as il x ~ 2A2/T = ilo/60, providing a 
very small contribution. 

Fig. [3] shows snapshots of the density profile for dif- 
ferent excitation times ranging from 13 to 17 ms. Ad- 
ditionally, we have calculated the mean angular momen- 
tum per atom, (L x ) — J dr^*L x '0, as a function of the 
excitation time. Using a = 1.6 and 7 = 0.02, our simula- 
tions show that (L x ) blows up after 15 ms of excitation. 
At this point, the condensate forms wavy patterns which 
develop to dark solitary waves which subsequently decay 
into several vortex pairs via the snake instability [29] . A 
more complex dynamics takes place after the first events 
of vortex formation, consisting in the generation of an 
undetermined number of vortices which characterizes the 
emergence of the turbulent regime. 

As a increases the nucleation of vortices occurs at ear- 
lier times, with a faster evolution to QT. This agrees 
well with the observations presented in the diagram of 
Fig. [5] It is important to mention that the time scale 
of the vortex events of the simulations, which is of the 
order of 10 ms, is consistent with the times observed in 
the experiment. 

These results demonstrate that the combination of ro- 
tation and translation is essential to produce vortices 
with lifetimes long enough to be observed after time-of- 
flight imaging. The instability on vortex generation may 
be associated with the Landau's argument for the criti- 
cal velocity [30] above which a fluid becomes turbulent. 
However, we must recall that, in general, the creation 
of vortices is not necessarily associated with the Landau 
critical velocity, as discussed in Ref. [3"T] . 

The dynamics of the BEC strongly depends on the 
strength of dissipation 7. If the dissipation is absent, 
no instability associated with the soliton creation occurs. 
Values of 7 between 0.015 and 0.025 are optimal for the 
generation of vortices and quantum turbulence. The sim- 
ulations cannot reproduce the whole observed diagram, 
since our system is a 3D gas. Nevertheless, good qualita- 
tive agreement with the experiment has been achieved. 



V. GRANULAR STATE 

When an overdose of energy is pumped into the turbu- 
lent cloud, the condensate enters in a new regime where 
all vorticity disappears and the system exhibits a new 
phase. In this case, the superfluid system demonstrates 
a kind of burst into grains. A typical image of this state 
can be seen in the insert (d) of Fig. [2J The resulting state 
is analogous to the Bose glass state, where the grains of 
condensate are surrounded by uncondensed atoms. The 
origin of this granulation can be understood as follows. 
The action of an external alternating field can be shown 
to be equivalent, on average, to the action of an external 
spatially random potential [32 . For an equilibrium sys- 
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FIG. 3: Snapshots of the BEC after different times of exci- 
tation. The left and the right columns show the 2D and 3D 
plots of the density profile, respectively. The colors range 
from red (high density) to blue (low density). 



tern, the granular condensate appears under the increas- 
ing amplitude of the external spatially random potential 
[261 133) . Similarly, for a nonequilibrium system, subject 
to the action of an alternating field, the granular state 
arises under the increasing value of the energy pumped 
into the system. The relation between the value of the 
oscillatory field amplitude, Vq, and the effective trap fre- 
quency wo is crucial for the occurrence of the correspond- 
ing states. This is connected with the relation between 
the characteristic trap length l n = l/^mwo, playing the 
role of the correlation length, and the localization length 
hoc ~ (wo/Vo) 2 ^Oj being an analog of the Larkin length 

E3I- 

When the energy pumped by the alternating poten- 
tial is small (Vo <C w ), there is an extended conden- 
sate filling the trap. This condensate can be strongly 
turbulent, but it still fills the whole trap. If the ampli- 
tude or modulation time of the alternating potential is 
increased, the condensate becomes strongly perturbed, 
eventually reaching the point where k oc ~ Iq. Under this 



condition, the condensate cannot be sustained as a whole 
and fragments into pieces. The nonequilibrium granular 
condensate can exist for the range u>q < Vq < Wo^/lo/o. 
When the pumped energy reaches the upper boundary, 
Vo ~ uio^/lo/a, then li oc ~ a and the granular conden- 
sate should be completely destroyed. This boundary was 
not reached in our experiment. Our observations are in 
agreement with the expected average size and number of 
the condensate droplets in the granular phase. 

The state, requiring very strong pumping and appear- 
ing after the granular condensate is destroyed, can corre- 
spond to nonequilibrium normal fluid in a chaotic regime 
|26j . This state can be similar to weak turbulence. The 
latter, as discussed in [35] . is not related to vorticity 
but is characterized by strong fluctuations of the den- 
sity distribution. The evolution from the granular to 
weak turbulent regime may be possible in the present 
system, which however, requires further experimental in- 
vestigations. Another interesting direction for experi- 
mental studies could be the process of relaxation of the 
nonequilibrium trapped system after the end of external 
perturbations. Finite quantum systems can demonstrate 
a number of peculiarities in the process of their equili- 
bration |3U1I37]. 



VI. CONCLUSIONS 

Turbulence presents one of the most challenging phe- 
nomena in physics. Even though this field has attracted 
high attention, its understanding is yet far from being 
perfect and a complete description of the problem is still 
missing. When turbulence occurs in superfluids it is 
known as Quantum Turbulence (QT) and presents many 
features which have no equivalent in its classical counter- 
part. In particular, the fact that the vortices are quan- 
tized makes the description of the turbulent flow simpler 
than in classical fluids. For this reason, it is believed that 
superfluids are excellent model systems for understand- 
ing such an important topic. To date, most of the work 
involving QT has been focused on the case of superfhiid 
Helium (He-4 and He-3). Nevertheless, the possibility of 
studying this phenomenon in a Bose-Einstein condensate 
(BEC) opens up new and exciting research opportunities 
mainly for three reasons which are not found in the case 
of Helium: 

(i) Condensates offer a very high degree of controlla- 
bility. In a BEC, it is possible to control the number 
of particles, interaction strength, and the confining po- 
tential where the atoms are trapped, making it an ideal 
system for generating QT in a highly controllable man- 
ner. 

(ii) Trapped atoms form weakly interacting systems, 
where the vortex cores are essentially larger than in su- 



perfluid Helium, making their observation much easier. 

(iii) The fact that the BEC is a finite-size system 
presents a new scenario for discovering new effects as- 
sociated with QT. 

In the present work we performed new studies of QT. 
In particular, we described the routes by which QT can 
be generated by means of an oscillatory excitation. We 
studied the evolution of a non-turbulent system to a tur- 
bulent one by varying the parameters of an oscillatory 
excitation applied to the sample. These observations, 
although peculiar for our system, are very important al- 
lowing for the general understanding of the conditions 
under which this phenomenon can be produced and in- 
vestigated. Besides the turbulent regime, when the high 
values of the parameters are employed, we also observe, 
for the first time, a different regime, which we have iden- 
tified as granulation of the BEC. In this case, the exci- 
tation is so high that the condensate breaks into small 
pieces surrounded by uncondensed atomic clouds. We 
summarized our observations in a diagram of excitation 
parameters, identifying the domains of different nonequi- 
librium states and clarifying the route by which these 
regimes could be reached. This diagram, being princi- 
pally novel, is peculiar of trapped atomic systems, whose 
size is finite and, therefore, their behavior can be dras- 
tically different from that of bulk superfluids, such as 
superfluid Helium. Together with our experimental ob- 
servations, we also provide a numerical model supporting 
our findings and allowing for a qualitative understanding 
of the evolution to quantum turbulence. We offered a 
plausible explanation of the arising granulated phase and 
suggested ways for its further investigation. 

The main results of the paper can be briefly summa- 
rized as follows. We have constructed a diagram showing 
the regions of different states realized under the com- 
bined variation of the time and amplitude of an oscil- 
latory excitation. The gradual evolution from a bended 
BEC to the regular vortex state, to the turbulent regime, 
and, finally, to granulated condensate is presented. Nu- 
merical simulations allow us to qualitatively explain the 
observations and to identify the requirements for real- 
izing this or that regime. The obtained diagram serves 
as a guide demonstrating the parameters that are nec- 
essary for experimentally producing different nontrivial 
nonequilibrium states of trapped atoms, such as turbu- 
lent condensates and granular condensates. 
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